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In recent years there has been an intensive search for Majorana fermion states in condensed matter 
systems. Predicted to be localized on cores of vortices in certain non-conventional superconductors, 
their presence is known to render the exchange statistics of bulk vortices non-Abelian. Here we 
study the equations governing the dynamics of phase solitons (fluxons) in a Josephson junction in 
a topological superconductor. We show that the fluxon will bind a localized zero energy Majorana 
mode and will consequently behave as a non-Abelian anyon. The low mass of the fiuxon, as well as 
its experimentally observed quantum mechanical wave-like nature, will make it a suitable candidate 
for vortex interferometry experiments demonstrating non-Abelian statistics. We suggest two exper- 
iments that may reveal the presence of the zero mode carried by the fiuxon. Specific experimental 
realizations will be discussed as well. 



I. INTRODUCTION 

Non-Abelian statistics^— has recently been the subject 
of intensive research driven both by its possibly profound 
impact on the field of quantum computation^— and by 
the search for its manifestations^. Among all mecha- 
nisms giving rise to such statistics, the route via spin- 
polarized p-wave superfluidity may be the simplest one. 
It was previously argued^r— that an Abrikosov vortex 
in a p-wave superfiuid can trap a zero energy Majorana 
fermion, being a self-conjugate "half" fermion. A pair 
of Majorana modes constitute a regular fermion, and the 
resulting non-local occupancies label a set of degener- 
ate ground states. Braiding of vortices results in mixing 
of these ground states, sometimes in a non-commutative 
fashion: it matters in which order multiple braidings are 
performed. The search for an explicit experimental sig- 
nal of the resulting vortex exchange statistics, as well as 
for the presence of Majorana modes on their cores, is 
currently on its way. 

In this paper we propose an experiment that probes 
Majorana fermions in Josephson vortices (fluxons). 
Josephson vortices are trapped in insulating regions be- 
tween super-conductors. For conventional superconduc- 
tors, they are described as solitonic solutions of the sine- 
Gordon equation moving with small inertial mass (esti- 
mated to be smaller than the electron mass) . In the case 
of topological superconductors, we find that such vortices 
bind a localized Majorana zero mode and would therefore 
behave as non-Abelian anyons, despite the fact that they 
lack a normal core. We show that the non-abelian nature 
of these vortices manifests itself in measurable transport 
properties of the Josephson junctions that house them. 

We start by showing that Josephson vortices traveling 
in a Josephson junction in a topological superconductor 
bind a single Majorana zero energy state (see Equation 
[BJ . We then discuss two experiments that can be used to 
measure the presence of these Majorana fermions. The 
first probes a thermodynamical property of a circular 
charge biased Jospephson junction (Figured]), by measur- 



ing the non-linear capacitance induced by the persistent 
motion of the vortex trapped in the junction^— . The 
second is an interference experiment of fluxons demon- 
strating Aharonov-Casher— oscillations (Figure [2]), sim- 
ilar in spirit to the one proposed ir*i£. 




FIG. 1. Aharonov-Casher effect in a long circular Josephson 
junction. The junction traps a single fluxon which is traveling 
around the ring propelled by a bias charge Q induced between 
the two ring-shaped superconductors. The energy spectrum 
of the junction is periodic in Q with periodicity e when <3> is 
increased to nucleate a vortex within the interior hole. Copper 
wires act as reservoirs of unpaired electrons. 



II. HAMILTONIAN OF A CIRCULAR 
JOSEPHSON JUNCTION 

We start by considering a circular Josephson junction, 
made of two concentric super-conducting annuli, sepa- 
rated by a thin insulator. We assume that the hole at the 
center of the inner super-conductor is of a size compara- 
ble to the superconducting coherence length, and encloses 
N v vortices. The Hamiltonian governing the junction 
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would be composed of three parts, H = H^ + H^ + H tun - 
The first, is related to the dynamics of the phase 
across the junction. For a Josephson junction of height 
h z , this part of the dynamics is derived from the following 
Hamiltonian (see e.gj 14 i 18 ) 



1 



1- -.2 1 



o +^2-(l-COS( 



(1) 



where cj> is the phase difference across the junction, n 
is the two dimensional (2D) density of Cooper pairs on, 
say, the inner plate, and a is the 2D density of the exter- 
nally induced charge. The first part of the Hamiltonian 
is the capacitive energy, the second the magnetic energy, 
and the third is the Josephson energy. The resulting 
equation of motion is the sine-Gordon equation, with the 



typical speed of light reduced to c 



d+2\ L _ 



(here d is 



the width of the insulating barrier and Al is the London 
penetration length). The Josephson vortex is a soliton 
described by this equation, with its typical size set by 
the Josephson penetration length Xj. We shall assume 
throughout that the circumference of the junction, L, is 
much larger than A,/, and that Al <C h z <C A,/. The 
parameter 1 = 16-r:j-^d(2X L + d)/h z . 

The second part of the Hamiltonian, H^, originates 
from the neutral protected edge modes of the topological 
superconductor, which give rise to its quantized thermal 
Hall conductance. The Hamiltonian governing these neu- 
tral modes is H^ = Hf + Hf with 
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where 1 and 2 refer to the two counter-propagating Ma- 
jorana edge modes. Here tpi(x) is a Majorana field, 
ipl(x) = ipii x )i x is the coordinate running along the 
Josephson junction, and v$ is the velocity along the edge. 
In terms of the electron's creation and annihilation op- 
erators c, (a;) , c\ (a;) and the phase fields of the two super- 
conductors on the two sides of the junctions <j>i{x) the 
Majorana fields may be expressed as 

^(x) = e^^^-^^c^x) + e- l<l} ^/ 2+ ™ x / L c\{x){?>) 

In this explicit form, the equation holds true for a spin- 
polarized p-wave superconductor. For a general topolog- 
ical superconductor, it is still the case that the Majorana 
field will acquire a minus sign going around the edge, 
in addition to one minus sign per each vortex enclosed 
in its path (i.e., one minus sign for every 2-7T winding 
of the phase fa). The boundary conditions of the fields 
ipi (x), ip2(x) therefore depend on N v , the number of vor- 
tices enclosed by the two super-conducting annuli. 

When the two superconducting islands are brought 
into close proximity, tunneling terms of the form 



c\(x)c2(x) translate to 

H tnn = 2im / dx ipi(x)ifa (x) cos(</>/2), (4) 



with 4>{x) = 4>2 {x) — 4>i (x) and m being a tunneling am- 
plitude. Writing Equations ([2]) and (HJ) compactly as a 
matrix equation, the Hamiltonian Hf + Hf + H tun be- 
comes 



H = J dx^ T {3 



iv^d x im cos (0/2) 
-im cos (0/2) —iv^d x 



(5) 



where ^ = (ipi,ip2) T is a spinor composed of the two 
counter propagating Majorana modes. The Hamiltonian 
possesses a symmetry under —¥ <j) + 2ti and — > a z ^ . 



III. BOUND MAJORANA MODE ON THE 
BACKGROUND OF A SOLITON 

A solitonic solution of Equation (fTJ), also known as 
a fluxon or a Josephson vortex, is a finite energy solu- 
tion which interpolates between two minima of the pe- 
riodic potential described by the Josephson term. For 
a long Josephson junction (L 3> Xj) it acquires the 

form 4> s (x) = 4arctan exp (^ x ^° ^j where xq is the po- 
sition of the soliton (see e.gi 14 i 18 ). In the following we 
solve Equation ([5]) in the background of a single soliton, 

explicitly plugging <j) s into <ft, and using cos ' > — 
- tanh ( ^ 



This would in turn result in a tunneling 

amplitude to cos whose sign is different on the two 

sides of the soliton. In light of the Jackiw-Rebbi mecha- 
nism, Equation ([5]) will now bind a zero energy mode at 
the position of the soliton Xq, 



7J = / dxf(x) [ipi(x) +sgn(m)V>2(aO] 



(6) 



In the limit of a long junction, L 3> Xj and L > »^/m, 
the shape of the Majorana mode is described by the lo- 

r / \ -| -|m|Aj/^ 

calized function f(x) — cosh I x x *° J , with 

N a normalization factor and xq the center of the Joseph- 
son vortex. The operator 7j is a localized Majorana 
fermion and the subscript J indicates that this mode 
is bound to a Josephson vortex. This mode satisfies 
~fj = j j. Indeed, the entire low energy spectrum of 
bound states can be extracted. Plugging the solitonic 
solution into Equation ([5]) and rotating the spinors ac- 
cording to ( \ , ] \[/ t he [famiHomaii can be 



V2 \ -1 
written conveniently as 



H = 



At 
A 



with 



A = —iv^dx — iW(x) 
A^ = — %Vif,d x + iW(x) 



(7) 
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The spectrum of this Hamiltonian possesses a supersym- 

metry with a superpotential W(x) — mtanh ^xj"^ • In 

particular, the zero mode is annihilated by either A or 
depending on the sign of the mass, while the other op- 
erator will not have a normalisable eigenfunction at zero 
energy. The rest of spectrum is doubly degenerate with 
A and A^ connecting the states of the doublet. Squar- 
ing the Hamiltonian in Equation (0) and linearizing the 
potential we get a shifted harmonic oscillator, and the 
spectrum of excitations above the zero energy mode is 



A, 



2n, with 



0,1,2, 



(8) 



Majorana fermions have to come in pairs, due to the 
properties of the BdG equations. The second Majorana 
fcrmion will be on one of the super-conductor edges, and 
will be denoted by j e , with e standing for edge. Its 
position will depend on the number of vortices N v in 
the inner hole. A zero mode is found on every edge 
that encloses an odd number of vortices. Hence, when 
N v is odd, 7 e will be localized on the edge separating 
the inner super-conductor from the vacuum. When N v 
is even, 7 e will be localized on the edge separating the 
outer super-conductor from the vacuum. The two Majo- 
rana fermions, j e and 7j, cannot shift from zero energy 
without hybridizing, and the super-conductor that sepa- 
rates them prevents that from happening. Note that in 
the absence of electron tunneling across the junction, the 
presence of these two Majorana modes is protected by an 
index theorem. Their wavefunction will be spread evenly 
around the two edges of the superconductor. When we 
next turn on electron tunneling, the Majorana mode 7j 
cannot disappear without hybridizing with the distant 
Majorana mode on the edge of the sample. Therefore it 
necessarily persists in the junction. Due to the effect of 
tunneling it gets localized around the center of the soliton 
as described by equation (J6j) . 

Due to the presence of the soliton in the junction, the 
boundary conditions are different for the two Majorana 
fields and depend on N v : ipi(x) = — (—l) Nv ipi(x + L) 
and tp2{x) = (— l) Nv ijj2(x + L). A rotation of the soli- 
ton around the junction shifts x by L and ^i(x) by 27r. 
These shifts have two effects. First, they multiply ipi by 
— (— 1)*^" and 1/J2 by (— 1)*^" due to the boundary condi- 
tions. Second, they multiply the off-diagonal term of ([5]) 
by —1 and hence multiply tpi(x) in Equation ([6]) by —1. 
Combining these two effects together, we see that 7j is 
multiplied by (— In view of this, the transformation 
7,7 — > (— 1) jv "7j may be understood as a consequence of 
the winding of one Majorana fermion around another for 
odd N v , and the absence of such winding for even N v . 
Remarkably, we find below that the capacitance of the 
junction in the two cases is different. 



IV. MANIFESTATION OF THE MAJORANA 
MODE IN THERMODYNAMICS AND 
TRANSPORT 

We now consider the energy of the long circular Joseph- 
son junction when it is biased by an external charge 
Q = 2eaLh z . For a "conventional" Josephson junction 
made of s-wave super-conductors this problem was stud- 
ied ini^. In that case the Hamiltonian consists of Equa- 
tion (TTJ) only. The magnetic and Josephson energies con- 
stitute the rest mass of the vortex. The charging energy 
constitutes its kinetic energy, and is the only component 
that depends on Q. This kinetic/charging energy may 
be understood in two ways. First, by viewing the vortex 



as a particle of mass M 



h 2 h z 



in a one dimensional 



ring, subjected to a vector potential with a set of 

energy eigenstates 
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E n {Q) = 



2ML 2 



Q 

2e 



(9) 



with n being the quantum number that quantifies the mo- 
mentum of the vortex. Second, by writing this energy as 

a capacitive energy of the form, E n (Q) = ^jj- — nj 

where now C is the effective capacitance (which for a 
short junction coincides with the geometric capacitance). 
The quantum number n is now identified as the number 
of Cooper-pairs charging the Josepshon junction, and the 
capacitance is C — M(2eL) 2 /(2nh) 2 . 

The spectrum of ((9|) is described by a set of parabolas, 
each characterized by a different integer value of n. When 
Q is increased so that the spectrum reaches a crossing 
point of two parabolas of n and n + 1, it becomes ener- 
getically favorable to tunnel a Cooper pair to decrease the 
charging energy, thus crossing to the next parabola. The 
matrix element required for this tunneling is generated 
by a small amount of disorder. The velocity of the mov- 
ing fluxon is proportional to the derivative of the ground 
state energy E gs (Q) = min ne z E n (Q) with respect to Q. 
In turn, the moving vortex generates a measurable volt- 
age via the Josephson relation between the inner ring and 
outer ring of the device, which will be oscillating as func- 
tion of Q with a periodicity of 2e. The periodicity of 2e 
is quite expected, given the energy cost needed to break 
a Cooper-pair. None of these results depend on the value 
of N v . 

The same spectrum shows a strikingly different behav- 
ior when the junction is made of two topological super- 
conductors. In particular, this difference manifests itself 
in the dependence of the spectrum on the charge Q and 
N v . The fermionic part of the Hamiltonian now has two 
degenerate eigenstates, which are two eigenstates of the 
operator ijj"f e , with eigenvalues ±1. When it comes to 
the spectrum of the Josepshon vortex, there is a crucial 
difference between the case where 7 e is localized at the 
edge of the outer super-conductor (even N v ) and that 
in which it is localized at the center of the inner super- 
conductor (odd N v ). As we saw before, in the latter case 



the winding of the soliton around the junction multiplies 
j j by — 1. It is easy to see that it multiplies 7 e by —1 
as well. These two operator transformations are imple- 
mented by the unitary transformation U — e la i^j~f e on 
the state of the system^. The phase a cannot be de- 
termined by this consideration, and is not important for 
what follows. The two ground states are eigenstates of 
U , and thus the application of U multiplies them by two 
phases, that differ by tt. The phase accumulated by the 
soliton encircling the ring affects its spectrum. When the 
energy (j^J) is viewed as the energy of a charged particle 
on a ring threaded by a magnetic flux, a tt phase shift 
corresponds to the introduction of half a quantum of flux. 
The spectrum (j^J is then changed to be 



Enj{Q) 



(2e) 2 
2C 
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where the / = ±1 correspond to the two eigenstates of 
U for a = 0. The charge Q' may be shifted with respect 
to the induced charge Q by a non-universal number that 
depends on a. The spectrum is now periodic with a pe- 
riod of e, rather than the period of 2e as was the case for 
the non-topological super-conductor, as well as for the 
topological super-conductor with an even N v . 

For the two sets of parabolas to be observed in an 
experiment when N v is odd, the system has to be able 
to switch between the two values of /, namely the two 
states ijjje = ±1. This switch requires the total number 
of electrons to change by one^ ' 19 ' 20 . To allow for that 
to happen, the junction should be weakly coupled to a 
reservoir of single electrons, i.e., to a metal. The ground 
state energy would be E gs (Q) = min neZj£Z2 E nJ (Q). 
The end result is that the voltage measured between the 
inner and outer edges of the superconductor would be 
periodic as function of Q (or Vq) with the periodicity 
now being e for odd N v and 2e for even N v . 

While this experiment probes a thermodynamical 
property, the next one is an interference experiment that 
probes electronic transport. The fluxons are generated 
near the bottom of the interferometer (Figure [5]), and 
driven by a super-current driven from left to right. The 
fluxon beam is split into two partial waves enclosing an 
island with externally imposed charge Q. Due to the 
Aharonov-Casher effect, the magnitude of the vortex cur- 
rent oscillates according to 



(11) 



where £ is the visibility of the oscillations and J„o the 
average current. By the Josephson relation, a measur- 
able voltage difference is created between the left and 
right ends of the superconductor, which would oscillate 
with changing Q. If the flux $ on the central island is 
increased so that another non-Abelian vortex is nucle- 
ated in the central hole, the interference pattern would 
be zero, £ = 0. These results are sensitive to effects of 
decoherence, which would be reflected in a reduction of 




FIG. 2. Vortex interferometry experiment based on the 
Aharonov-Casher effect adapted to Josephson vortices. A 
superconducting wire creates a circulating magnetic field act- 
ing as a source for the entrance of Josephson vortices into the 
sample. An applied super-current drives the vortex along one 
of two paths circumventing an island towards the top of the 
sample. A charge Q enclosed in the island controls the inter- 
ference term via the Aharonov-Casher effect. When the flux 
$ nucleates a vortex in the central region, the interference 
term would be obliterated. 



the value of £ as the temperature is increased. How- 
ever, the fluxons are quite light, and it was predicted 21 
and experimentally demonstrate d 15 ' 22 that they can ex- 
hibit quantum behavior 23 ! 24 . Also, the gap to the next 
fcrmionic state on the background of the soliton has dif- 
ferent parametric dependence as compared to Abrikosov 
vortices. 

For a soliton in a topological superconductor there are 
two main energy scales that control decoherence effects. 
The first is related to the bosonic part H^, while the 
second stems from the fermionic part . 

The spectrum of the bosonic Hamiltonian contains 
non-topological excitations with which the fluxon can in- 
teract. These are quantized phase oscillations, or plas- 
mons. Written in terms of collective coordinates for the 
fluxon, the Hamiltonian would be similar to the Caldeira- 
Leggett Hamiltonian of a quantum particle interacting 
with a bath. However, unlike the Caldeira-Leggett mech- 
anism, here the plasmon spectrum is gapped, with the 
plasmon gap given by hu> p — hc/Xj. At energy scales 
below this gap (estimated to be a few Kelvin) the cou- 
pling to plasmons is exponentially suppressed, and the 
internal dephasing thus mostly avoided^ 5 - over lengths 
much larger than the Josephson penetration length. This 
lies at the origin of the prediction of quantum behav- 
ior of Josephson vortices in long Josephson junctions. 
In particular, scattering of the fluxon on static inho- 
mogeneities in the junction would be mostly elastic due 
to the presence of this plasmon gap 21 . Tunneling of a 
fluxon through barriers was predicted to be enhanced 
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by the presence of the plasmons^i and experimentally 
observed 2 ^. 

The second energy scale is unique to topological su- 
perconductors, and is related to the presence of higher 
energy fermionic states above the zero energy mode car- 
ried by the soliton, as described by if ^ and the fermionic 
tunneling terms, see Equation ([5]). The presence of 
these states would result in decoherence of the non- 
Abelian mechanism. The gap to the first state would 
be hui n — y/2hv^m/Xj. To avoid these effects, the tem- 
perature should be lower than this gap as well. 



V. PROPOSED REALIZATIONS 

There are several systems in which the ideas presented 
above may be implemented in practice. The most rele- 
vant for our purpose is a topological insulator whose sur- 
face states become superconducting by proximity to an 
s-wave superconductor, and are driven into an effective 
p-wave pairing for a single fermionic species. The second 
is a hybrid structure of three materials: a semiconduc- 
tor with a magnetic material and a superconductor layer 
placed on top of it, as discussed in^ 6 -. The third is the 
perovskite material Sr2Ru04 (SRO) which becomes su- 
perconducting at temperatures below 1.5K, and for which 
there is growing evidence that it realizes an unconven- 
tional pairing of a p x ± ip y (spin-triplet) form. 

The surface state of a topological insulator is described 
by the Hamiltonian ffo = Vpip'z ■ cr x (— iV)?/> (where 
-0 = (ipf, ^i) T is a spinor composed of the electronic op- 
erators associated with the spin components) . The chem- 
ical potential is tuned to the Dirac point. Superconduc- 
tivity is induced on the surface state by the proximity 
effect^, represented by an extra term A^l^t + h.c. to 
the Hamiltonian. The superconducting layer of height 
h z is deposited as described in Figure [U At the thin in- 
sulating layer between the two annuli, as well as in the 
inner and outer holes, an insulating magnetic material 
of constant magnetization should be deposited, resulting 
in a Zeeman term of the form Mip^a z tfj in the surface 
Hamiltonian. This breaks the time reversal symmetry at 
the edges and by that chooses a single chirality for the 
flow of the neutral chiral edge modes. The dynamics of 
the soliton will be largely determined by the s-wave su- 
perconducting layer, while a zero energy Majorana mode 
will be trapped by the soliton on the surface state of the 
topological insulator. Other than these material specific 
details, the rest of the arguments in the paper can be 
applied with no further changes. 

To estimate the experimental parameters we take typ- 
ical values Xl = 0.1/im and Xj = 30/im for an s-wave 
superconducting layer with h z = 0.5/im and an insu- 
lating region of width d — 20A. For these parameters, 
the dimensionless parameter 1 ~ 0.01, and the veloc- 
ity of light is c = 0.1c. We take the neutral edge ve- 
locity to be tty = 10 5 m/s2£ (which coincides with the 
surface state Fermi velocity close to the Dirac point), 



and the tunneling between Majorana edge states to be 
0.025meV. The plasmon gap is E p ~ 5K2£, while the 
intra-core gap is given by direct substitution to equation 
([S}, E n ~ 120mK. The junction charging energy is E c = 
(2e) 2 /2C ~ 250(Aj/L) 2 mK. In order to observe quan- 
tum phenomena, we need the temperature to be smaller 
than all these energy scales, T < min {E n , E c , E p }. The 
operating temperatures are therefore in the range of 
10- 100 mK. 

The implementation of the Josephson junction in a 
hybrid structure of semiconductor quantum well cou- 
pled to an s-wave superconductor and a ferromagnetic 
insulator 26 is similarly straightforward. The supercon- 
ducting layer is deposited as in Figure [TJ while the fer- 
romagnetic insulator layer will be deposited throughout 
with no restrictions, breaking the time-reversal symme- 
try for the quantum well. Another option is to use a 
quantum well with both Rashba and Dresselhaus spin- 
orbit coupling, and deposit the superconducting layer as 
before, with an external magnetic field applied in the in- 
plane direction^. Both cases will result in an SNS type 
Josephson junction. 

The implementation using SRO is different in several 
respects. First, SRO is not spin polarized. As a con- 
sequence, in bulk SRO Majorana modes are carried by 
half-vortices, that are made of a ir rotation of the pair- 
ing d-vector glued together with a tt phase winding of 
the order parameter. These half-vortices were recently 
observed 2 ^ in mesoscopic samples of the type that may 
be useful for the present context. Similarly, Majorana 
modes in a long Josephson junction would require half- 
fluxons. Second, SRO is a three dimensional material, 
made of two dimensional layers. The physics described 
above for the Majorana modes carried by semi-fluxons 
decouples into different layers, one Majorana edge per 
layer. The multitude of Majorana modes does not af- 
fect the even N v case. In the case of odd N v , the factor 
7e7,7 in U is replaced by Y[f=i le,ilJ,h with the index i 
numbering the layers, and L z being the number of layers. 
The unitary transformation still has two eigenvalues that 
differ by a minus sign, and a transfer of a single electron 
facilitates a transition between the two states. Thus, as 
long as the zero energy states do not split, the multi-layer 
case will not be different from the single layer one. 

Tunnel coupling between the layers may split the Ma- 
jorana zero energy modes. However, a tunneling term 
between the i,j layers is of the form 

ifcoupiing = 2it J dxipi(x)ipj(x) sin [(fa - (j>j)/2) ,(12) 

with z, j being the layers involved in the tunneling. This 
tunneling is suppressed when the phase difference be- 
tween the layers vanishes. Fluctuations of 4> are massive 
due to the presence of the inter- layer Josephson coupling, 
and thus we may expect the Majorana modes not to split 
(see also^ -). 

Third, SRO may take two possible forms of p-wave 
pairing, commonly denoted by p x ± ip y . The analysis 
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above assumes that both super-conductors are of the 
same pairing form. If the converse is true, the shape 
of the Majorana zero mode is different, but the rest of 
the analysis is unaffected. We analyze this case in the 
supplementary. 

We note that SRO suffers from the presence of a small 
mini-gap ~ A 2 /Ep for quasi-particle excitations in the 
core of bulk vortices and in edge states. This sets a con- 
straint of about lOmK on the maximum temperature. 
For the case of a topological insulator tuned to the Dirac 
point this problem is avoided. 

VI. SUMMARY 

In summary, we predict that phase solitons in a long 
Josephson junction embedded in a topological supercon- 
ductor would carry a localized Majorana zero mode. 
Consequently, these solitons would constitute anyons 
with non-Abelian exchange statistics. Exploiting the 
quantum nature of these solitons we suggest two ex- 
periments that can reveal the presence of the Majorana 
modes. One experiment involves voltage oscillations of a 
topological superconducting capacitor, realized by a cir- 
cular Josephson junction hosting a single fluxon. The 
other experiment involves interference effects of a fluxon 
beam. 
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VII. SUPPLEMENTAL INFORMATION 

A. Electron tunneling between Majorana edge 
states 

In this section we derive the low energy form of 
the electron tunneling term across a Josephson junction 
made of two topological superconductors separated by 
an insulator. We take the topological super-conductors 
to be p x ± ip y spin-polarized superconductors in the weak 
pairing phase, where p, > 0^2. We take the insulator to be 
a region where /i < 0, and find it convenient to allow for 
a pairing term also in the insulating region. We do not 
expect the physics to be altered much by the small oc- 
cupations induced by the pairing term in the conduction 
band of the insulator. Also, we do not expect the results 



to depend on the microscopic details - they should apply 
for any topological superconductor in the same univer- 
sality class. As we describe below, in the absence of elec- 
tron tunneling across the insulator the model predicts 
the presence of two spatially separated Majorana edge 
modes with a tunnel barrier between them. The modes 
are counter-propagating for a junction separating two su- 
perconductors of the p x + ip y type, and co-propagating 
for a junction separating &p x + ip y super-conductor from 
a Px — ipy one. 

We start by reviewing the case of a single straight in- 
finite edge state, worked out in details i n 10 i 31 ' 32 . The 
Hamiltonian density of a p x + ip y superconductor in the 
presence of an edge going along the x-axis (y — 0) is 
C^HC/2, where C — (c, c^) T , c is the electron operator, 
and 

»-(3_£ir). m 

with 

H a = -±-V 2 + V{y)-ti, (14) 

A = iv A e l4 > (r)/2 (d x - id y y^ r)/2 . (15) 

Note the following symmetry of the Hamiltonian, 
a x H*a x = —H, having the consequence that for any 
solution xe — ( u e, v e) T with energy E there exists a 
solution X-e = &xXe with energy — E. Define jl(y) = 
fj, — V(y) the electro-chemical potential, which we take to 
be negative on one side of the edge (the insulator), and 
positive on the other (the topological super-conductor). 
Linearizing the Hamitonian near the edge, we get for 
H> = e --J <T ^( r )_H e 5 <T ^( r ) ) 

H' = -<j z {p{y) - v A a x d y ) + iv A a x d x . (16) 

The Hamiltonian factorizes into x and y components. 
The part that depends on y is an effective one- 
dimensional Dirac Hamiltonian with a mass term that 
changes sign at the edge, and therefore contains a zero 
mode. We assume that the potential is sharp enough so 
that the higher energy states can be disregarded. We can 
now write the Nambu spinor projected into this zero en- 
ergy state, C e (r) (with the subscript standing for "edge") 
as 

C e (r) = e ^«e^ # Ky ' )dy ' ( ^ \ ^(x), (17) 

where the effective Hamiltonian for the Hermitian field 

ip(x) is 

H" = ±iv A d x , (18) 

and ± is the eigenvalue of a x . As is apparent from equa- 
tion (fT6|) , the eigenvalue of a x is uniquely chosen by the 
behavior of the electro-chemical potential fi at y — ¥ ±oo 
to make the zero-mode in the y-direction normalisable. 
This effectively locks the direction of the propagation of 
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the edge to the spatial profile of sgn(/i(y)). The operator 
ip{x) can now be expanded as 



(19) 



where ipk satisfies ipt — ip-k an d L is the length of the 
junction. The energy is related to k via E = ±v&k, 
identifying v^, — «a- 

It will prove useful to generalize the result to a straight 
edge forming an angle 9 with the z-axi o 31 i 32 , for two rea- 
sons. First, it gives a handle on the direction of propaga- 
tion of the edge, including inverting its chirality. Second, 
when we consider a circular edge, we can locally describe 
the edge as being flat, with the parameter 9 affecting the 
boundary conditions through an explicit monodromy. To 
proceed, we denote by X the direction along the edge, so 
that (x, y) = X(cos9, sin 9), and define Y as the perpen- 
dicular direction, 



C e (R) 



e 2 



1 

±1 



(20) 



where R — (X, Y). Note that a rotation by an angle 9 = 
7T reverses the chirality of the edge, connecting the two 
orthogonal chirality states represented by eigenvectors of 
a x . Also, for a large circular edge we can treat the term 
which depends on 9 as a rotating "frame" , leading to anti- 
periodic boundary conditions. Therefore, when an even 
number of vortices is enclosed by the edge k is quantized 
to be k = ^jj{n + 1/2) (n g Z), while for an odd number 
of vortices k — ^-n , so that the spectrum contains a 
zero modei^. 

Suppose we have two edges, at y = ±a/2, where p,(y) < 
for —a/2 < y < a/2 and ft{y) > otherwise. Assume 
the form of the order parameter is p x + ip y throughout 
the sample. We can write the Nambu spinor projected 
on the two edge states as 



C e i = e 



4 K /a Hiv'W e 2 



i$2 

e 2 



(21) 



Hence, electron tunneling between edges is of the form 

0i (x) - 4> 2 {x) 



im / dxipi(x) ip2(x) cos 



(22) 



In the case that the insulator separates a p x + ip y su- 
perconductor from a p x — ip y superconductor, for the 
Px ~ iPy sector there will be an extra minus sign preced- 
ing the term containing d y in equation (fTo| . This would 
flip again the eigenvalue of a x with respect to the orien- 
tation, so that C ej i will be redefined to be 



Ce,l 



x Ia/2Mv') d v' 



e 2 



M*), (23) 



leading to the following form for the tunneling term 
im f dx ip i (x)->Jj 2 (x) sin * 



(24) 



B. Solutions for the Majorana mode for different 
cases 

In this section we find zero energy solutions of the neu- 
tral sector of our main Hamiltonian. Written in a general 
form, the Hamiltonian is 



H = 



ivipd x iW±(x) 
-iW±(x) ±iv^,d x 



(25) 



where W+(x) — msin(0/2), W~(x) = m cos(</> / '2) , and 
± denote the relative chirality of the two edge states. To 
solve for zero modes, we write H\ = where \ — (f,9) T - 
Squaring the Hamiltonian, we get 



g2 _ dsW±(x) ft // „..,,, 



W± (x) m z 



g(x) = (26) 



Here rj = ^ with £ a length scale related to the 

potential W±(x), and x — xjl. A general solution 
for co-propagating Majorana states is (define $±(i) = 
If W ± {x')dx>) 

f+ (x) = - A sin [$+ (x)] + B cos [$+ [x)] , 

g+(x) = A cos [$+(5)] + Bsin[$ + (i)] . (27) 

For counter-propagating Majorana states one gets 

/_ (J) = -A sinh [$_ (£)] - B cosh [$_ (x)] , 

g- (x) = A cosh (x)] + B sinh [$_ (x)] . (28) 

The resulting Majorana fermion acquires the form, 

7J,± = [dx [f± (x)ih (x) + g± {x)i> 2 {x)\ . (29) 



In the body of the paper we derive the zero energy 
modes for the case of a long Josephson junction separat- 
ing two p x + ip y super-conductors, where the two Ma- 
jorana modes propagate in opposite directions. In the 
following, we shall write down explicitly the solutions for 
two other case, that of co-moving Majorana modes, and 
that of a short Josephson junction. 



Co-propagating Majorana 



states 



The Hamiltonian of the neutral modes for a p x + 
iPy I Px — ipy j unction is different from a p x + ip y / p x + ip y 
junction for two reasons: First, the Majorana edge states 
are propagating in the same direction, and second, the 
cos(</>/2) in the mass term is replaced with sin(0/2), with 
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p + ip 



p + ip 




FIG. 3. A Px+ipy /p x +ipy junction versus a p x +ip y / p x —ip y 
junction. The direction of the propagation of the Majorana 
edge states is marked by arrows. Dotted lines indicate elec- 
tron tunneling. 



4> the phase difference between the two superconductors. 
The resulting Hamiltonian in the presence of a soliton is 



H = 



ivtj)d x 
-imsech (^jj^J 



im sech 

iv-,pd x 



(30) 



For a closed long circular junction the zero energy solu- 
tion to that equation is of the class presented in equation 
(g7J) , witn 



$+(x) = 27/arctan 



tanh 



2Aj 



(31) 



where r\ = m\j /v^ and — ^ < x < 4 . The two functions 
/+ and g+ should also satisfy the boundary conditions 
imposed by the magnetic flux configuration. For an odd 
number of vortices within the central hole / is periodic 
while g should be anti-periodic, hence A — 0. Similarly, 
for an even number of vortices B = (In the absence of 
a fluxon both / and g have the same boundary condi- 
tion, hence both A and B must vanish, and there is no 
zero energy solution). In contrast to the case of counter- 
propagating edge modes, where the Majorana state is 
localized close to the center of the soliton, in the case of 
co-propagating edge modes it is delocalized throughout 
the junction. However, the profile of the fluxon imposes 
a sign change on the part of the Majorana mode's wave- 
function which is bound to the anti-periodic edge. Hence, 
when the fluxon is adiabatically carried around the junc- 
tion, the boundary conditions result in the usual trans- 
formation applied to this delocalized Majorana mode, i.e. 
U -+ (-1)^7.7- 



limit A j ^> L (a short Josephson junction). In this limit, 
the phase interpolates linearly between and 27r, so we 
can write <fi = + n (taking the center of the soliton 
to be at x = 0). For a junction of counter-propagating 
edge modes, the neutral sector of the Hamiltonian takes 
the approximate form 



H 



iv^d x — zmsin(^) 
imsin(^7) —iv^d x 



(32) 



where r/ = mL/v^p. The solution is of the class presented 
in equation (|28[) . with 



*-(*) = " cos ( T ), 



(33) 



where we have dropped a constant term by redefining 
A and B in equation (|28|) . Most of the weight of the 
wave-function is therefore present in the vicinity of the 
soliton. 

For a short junction where the edge modes are co- 
propagating the solution is of the class presented in equa- 
tion (|2T|). with 



(34) 



For this case, the Majorana solution is uniformly spread 
throughout the junction, with its weight shifting between 
the inner and outer edges to conform with the boundary 
conditions. Here it is straightforward to see explicitly 
the effects of the boundary conditions. When x is in- 
creased by L, we get that §±(x + L) = — <&±(cc) and, 
consequently, the odd functions (sin, sinh) in equations 
(J2"7 |) - (|2"5)) change sign, while the even functions (cos, cosh) 
stay the same. Due to the boundary conditions imposed 
on / and g, one periodic while the other anti-periodic, ei- 
ther A = or B = 0. When combined with <f>i — > 4>i + 2tt 
(originating from the effect of the fluxon on the phase of 
the superconductor enclosed by its path, and leading to 
ipi — > —ipi), this would result in the usual transformation 
rule for the Majorana operator defined by this solution, 
i.e. 7j -> (-l) Ar "7j. 

In summary, independent of the size of the soliton rela- 
tive to the size of the junction and of the relative direction 
of motion of the two edge modes within the junction, a 
winding of the soliton around the junction results in the 
transformation 7,7 — > (— 1) "7,7. Thus, at least to the ex- 
tent that we deal with a single fluxon in the junction, the 
effects associated with the Majorana zero energy state in- 
duced by the fluxon in the junction are independent of 
these details. 



C. Classical fluxon propagated by a super-current 



2. The case of a large soliton 

In the main text we consider the case that Aj <C L 
and — < Aj. Here we will discuss the solution in the 



In this section we consider the bound states on a clas- 
sical fluxon moving at velocity v^. We assume that the 
junction is oriented along the x-direction. To move the 
fluxon, we apply a super-current along the y-direction. 
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For low values of the applied current / the attained ve- 
locity of the fluxon is linear in / with a coefficient that de- 
pends on the microscopic details of the junction^. When 
the fluxon moves at a velocity v^, it generates a voltage 
along the perpendicular direction given by 



V = (pov^/L, 



(35) 



where L is the length of the junction and <fio = h/2e is the 
superconductor flux quantum. This behavior generates 
the typical I — V characteristics of the junction. 

The Hamiltonian for the neutral states in the back- 
ground of a moving soliton is given by 



iv^d x iW±{x-v<j ) t) 
-iW± {x - v<f,t) ±iv^d x 



*±,(36) 



where if? = (-01, ip2) T ■ We solve this equation by going to 
the frame which moves with the soliton. Defining x' = 



x — v^t to be our new coordinate, 
d t = d t > - v$d x 



d — d i 



we consequently get 



iv^dx' iW±(x') 
-iW±(x') ±iv^d x ' 



(37) 



*' ± (38) 



For co-propagating Majorana edge states the solution 
is trivial, since ity is simply replaced by + v$. The 
more interesting case is when the Majorana edge states 
are counter-propagating, and the Majorana fermion is 
localized. Then, when solving for the zero mode, we 
can square the Hamiltonian and get the result appearing 
in the main text for the form of the zero energy wave- 
function, but with the parameter r\ = rnXj/v^ replaced 

by r]/\Jl — (v^/v^p) 2 . The size of the Majorana mode 
wave-function therefore decreases with increasing speed 
of the soliton, with a Lorentz-like dependence on its ve- 
locity and an effective speed of light which is set to the 
Majorana edge velocity v^p (the soliton decreases in size 
as well, with the effective speed of light being c >■ tty). 
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